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The prediction of manufacturing yield of an electric circuit, given 
the tolerances and statistics of the components, is a straightforward 
procedure using Monte Carlo tolerance analysis. The inverse problem, 
namely, specifying the component tolerances that produce the cheap- 
est network, is more difficult and has received little attention in the 
literature. In this paper an algorithm is presented that solves this 
problem for a significant class of networks. 

The algorithm is limited to circuits for which one hundred percent 
yield is sought and whose components are statistically independent, 
i.e., discrete circuits. The one hundred percent yield assumption is used 
to reduce the number of possible tolerance choices. A variation of the 
branch and bound strategy that is shown to be particularly efficient 
is then used to make an optimum selection. Two-dimensional per- 
formance contours, worst-case, and Monte Carlo computations are 
also used as part of the procedure. 

An example circuit is demonstrated for which the algorithm pro- 
duced a tolerance assignment substantially cheaper than that pro- 
jected by the designer. 

I. INTRODUCTION 

Monte Carlo tolerance analysis has proven to be a useful tool in 
evaluating the effects of component tolerances and environmental 
variations on electrical circuit performance. The method involves 
"constructing" samples of the circuit inside the computer using ele- 
ment values that obey the manufacturing statistics, analyzing these 
samples, and forming empirical distributions of performance. One 
common outcome of the process is the prediction of yield. 

Monte Carlo tolerance analysis, henceforth referred to as TAP 
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(Tolerance Analysis Procedure), is an open-loop structure. If we ex- 
amine the way in which it is used, we find that for discrete circuits 
the designer supplies the TAP program with a set of component toler- 
ances he has chosen based on breadboard measurements, linear sensi- 
tivity or worst-case analysis, or some other approximate technique. 
At the conclusion of the TAP run, he observes yield and is faced with 
one of two situations: 

(i) Yield is too low. With this result the designer knows he must 
change his tolerances. Unfortunately, TAP gives him little 
information as to which tolerances to change and by how much. 

(u) Yield is adequate. Here the designer may be satisfied by the 
design but he obtains little help in determining whether a 
cheaper (looser) set of tolerances might not give equally satis- 
factory yield. 

TAP, then, is open-loop in that it is a tool for predicting yield given 
a set of tolerances ; both the initial set of tolerances or any changes to 
that set must be provided by the designer without assistance from 
TAP. (It is interesting to note that for integrated circuits the designer 
has little freedom in specifying tolerances or tracking. In this case 
TAP is used to check that designs themselves are adequate, given the 
component tracking that can be achieved in integrated production.) 

This paper discusses an algorithm for closing the TAP loop. The 
program embodying this algorithm will accept a circuit topology, nom- 
inal element values, a figure of merit for circuit performance and data 
relating element cost to element tolerance, and specify the set of 
tolerances that will produce the cheapest manufactured network that 
has 100 percent a priori yield. (By a priori yield we mean the per- 
centage of those manufactured circuits all of whose components are 
within tolerance limits and which are wired correctly that meet specifi- 
cations. 100 percent a priori yield implies that any manufactured 
circuit that fails to meet specifications has either a component outside 
of tolerance or a wiring error.) The techniques described are appli- 
cable at present only to discrete circuits, i.e., circuits where the 
parameters are statistically independent. Work is in progress to treat 
circuits with correlated parameters. 

The method described below is a combination of two-dimensional 
performance contours, 1 a method the author has since learned is a 
variation on the branch and bound technique, 2 quasi-worst-case com- 
putations and repetitive TAP computations. 
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II. THE CHEAPEST NETWORK 

In this paper the cheapest network is defined as the network with 
100 percent a priori yield such that the sum of the costs of the com- 
ponents comprising that network is not greater than the sum of costs 
for every other network with 100 percent a priori yield. We are con- 
cerned with a network of fixed topology and fixed nominal element 
values for which we want to choose a tolerance vector, r, each of whose 
elements, n , is the tolerance chosen for circuit parameter, p t . Associ- 
ated with each parameter tolerance, Tf , is a cost, CiM, where Cd-) is 
the cost function of the ith circuit parameter. 

For each parameter the designer has a limited number of discrete 
tolerances from which to choose. Therefore, we will consider a finite 
universe of possible realizations, each realization being characterized 
by a unique tolerance vector. We will number the realizations in the 
universe 1, 2, •••, I, and refer to each realization by its associated 
tolerance vector t*. We next define the set, 

a = {k | realization k has 100 percent a priori yield} 

and the associated set, 

a T = {r* | k t a}. 

Then a cheapest network has tolerance vector t* satisfying the condi- 
tions: 

(i) r*ea T , 

iii) E C,(r*) =g £ C,(t*) ; for all k t a, 

where the circuit under study has n components. 

One could imagine other criteria by which the cheapest network 
might be defined. In particular, one might find a network cost function 
that traded off yield against cost of repairing, component salvaging, 
or discarding networks. However, consideration of such tradeoffs 
raises serious questions about manufacturing practices, not to mention 
the question of allowing sub-marginal circuits to be installed in the 
field. Because of these questions, the following discussion is limited to 
circuits for which 100 percent a priori yield is required. As we shall 
see, this assumption is central to the proposed tolerance specification 
method. 
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III. REGIONS OF ACCEPTABILITY AND POSSIBILITY 

It is assumed that the circuits with which we are dealing have a 
scalar performance, J, that is a computable function of the circuit 
parameter values; a nominal design with performance measure, Jo \ 
and an allowable degradation from nominal, e, such that (J + e) 
marks a boundary separating acceptable circuits from unacceptable 
ones.* 

We now examine the n-dimensional parameter deviation space, i.e., 
a space whose origin is defined as the nominal parameter vector and 
each of whose axes is the percent deviation of one parameter from 
its nominal value. We postulate the existence of a connected region, 
including the origin, such that all circuits built with parameter values 
represented by points inside the region are acceptable and those out- 
side the region are unacceptable. We call this the Region of Accepta- 
bility, jR.4 - 1 We now notice that in this same space, each tolerance 
vector, r k , is associated with an n-dimensional parallelepiped that we 
call the Region of Possibility. 1 In other words, all circuits built of 
components with tolerances specified by r k have parameter deviation 
vectors inside the parallelepiped we are calling the Region of Possi- 
bility, Rp . 

To clarify these ideas we extend an example presented in Ref. 1. 
Figure la shows a simple voltage divider. Its transfer function is given 
by T = 1/(1 + R x /R 2 ), and its input resistance is R = R t + R 2 . With 
nominal one-ohm resistors, the nominal transfer function and input 
resistance are 0.5 and 2.0 respectively. Suppose the design specifications 
call for 0.46 ^ T ^ 0.53 and 1.85 £ B ^ 2.15. Then the Region of 
Acceptability (R A ) shown in Fig. lb is bounded by four lines, each of 
which is the locus of all points producing networks exactly satisfying 
one of the four specification limits. Figure lb also shows the Region 
of Possibility (R P ) when the voltage divider is built with 5 percent 
components. Notice that R P intersects the parameter deviation axes 
at ±5 percent. Notice further that all possible networks are acceptable, 
i.e., R P is totally contained within R A . In Fig. lc, a new R P that is the 
result of using a 10 percent resistor for Ri and a 3 percent resistor 



* See Ref. 1 for examples of transcribing classical loss, phase, frequency specifi- 
cations into /, e representation. 

t Here we are assuming that all components are 100 percent tested, i.e., a batch 
of 15 percent components has all components within ±15 percent of nominal 
value. We are further assuming that tolerances place symmetric limits on com- 
ponent values. The asymmetric case can be handled but it makes the argument 
obtuse. Lastly, we are assuming all components are statistically independent. 
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Fig. la — Voltage divider example. 

Fig. lb — Regions of acceptability and possibility for voltage divider (5 percent 
components). 

Fig. lc — Regions of acceptability and possibility for voltage divider (Ri a 10 
percent resistor, R 2 a 3 percent resistor). 



for R 2 is superimposed over R A . Notice that some networks in the 
second quadrant fail, i.e., they are inside R P but outside R A . Hence, 
the tolerance vector (10, 3) is not a member of a, and cannot be used 
for the example network. 

The conclusion we may now reach is that in order for a network to 
have 100 percent a priori yield we must have R P Q R A . Further, 
because we are restricting our attention to networks that have 100 per- 
cent a priori yield, our task is to choose the cheapest tolerance vector 
that will satisfy the condition, R P C R A . 

If we could characterize the Region of Acceptability, the task of 
choosing a tolerance vector that produces the cheapest network would 
be reduced to the problem of finding the "largest" parallelepiped (R P ) 
that is contained within R A . (Notice that our assumptions of 100 per- 
cent a priori yield and 100 percent component testing obviate the need 
to consider component distributions. All we are concerned with are the 
limits of component deviation.) Unfortunately, characterization of R A is 
an impractical task; we are trapped by "the curse of dimensionality." 
Even if R A were as simple a region as a hypercube, it would, for a net- 
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work of say 15 parameters, have over 32,000 vertices. The number of 
function evaluations required to locate these vertices would be exceed- 
ingly large. When one considers that the (n — 1) dimensional surfaces 
bounding R A are likely to be quite complex, the magnitude of the 
problem becomes staggering. 

The algorithm that is used to select the tolerance vector that will 
produce the cheapest network, r*, begins with a consideration of 
two-dimensional subspaces of R A .' We will require that two-dimensional 
subspaces of R P be inside the corresponding subspaces of R A . We will 
then use a series of techniques to insure that R P Q R A . 

IV. FINDING THE CHEAPEST NETWORK 

The algorithm for finding a cheapest network is broadly outlined 
in Fig. 2. Below we take up each of the major blocks in detail. 

4.1 Pairwise Constraints 

The first step in obtaining pairwise constraints on component tol- 
erances is to compute the minimum intercept 1 for each parameter. 
Call this Li for the zth parameter. Briefly, this is the minimum devia- 
tion each parameter can undergo, keeping all other parameters fixed 
at nominal value, before the network specifications are violated. 

Next we consider the (J) performance contours 1 for the n parameter 
network. These contours are the boundaries of the two-dimensional 
subspaces of R A , holding (n — 2) parameters fixed at nominal. Each 
contour determines what tolerance pairs are allowable for the variable 
parameters. In short, each tolerance pair defines a rectangle. The 
pair is allowable if that rectangle is within the performance contour. 

The designer is presumed to have specified a discrete set of obtainable 
tolerances for each component. Call this set S t for the ith component. 
We will number the members of *S, in descending order and refer to 
them as r iX , t, 2 , • • • , r im , where m, is the number of obtainable 
tolerances for the ith component. We will form a table for each pair 
of parameters, say (r, s), of the form: 



Pr 


P> 


Tr.i 


T,, a 


Tr.m r 


T..y 



For ease of subsequent computation, the parameter on the left will be 
the one with the smaller minimum intercept. Not all m r entries neces- 
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YES 



FIND PAIRWISE CONSTRAINTS 
IMPOSED BY Ra 2 FOR 
ALL TOLERANCE PAIRS 



INITIALIZE: CHEAPEST 
NETWORK TO ZERO COST 



FIND NEXT CHEAPEST 
NETWORK SATISFYING 
PAIRWISE CONSTRAINTS 



CALCULATE "WORST CASE" 
PERFORMANCE 



DOES 

"WORST- CASE 1 
SATISFY 

CONSTRAINTS ? 



YES 



RUN TAP WITH WEIGHTED 
DISTRIBUTIONS UNTIL A 

FAILURE OCCURS OR 
UNTIL 300 HAVE PASSED 



HAS A FAILURE 
OCCURED ? 



PRINT OPTIMAL 

TOLERANCE VECTOR 

-* 



Fig. 2 — Tolerance specification algorithm. 

sarily exist in the table. First we eliminate those tolerances that are 
larger than the minimum intercept, since they will clearly produce rec- 
tangles not contained within the (r,s) contour (r r .„). Then for the 
first T r ,j ^ \L r \, we find the largest element of S s that produces a rec- 
tangle inside r r „, , and this tolerance is entered into the table opposite 
T r j . If a line of the table reads 



Tr.i 



T'.k 



and if parameter r has tolerance r r ., , parameter s may have r., k or 
any smaller tolerance. We proceed through the elements of S r finding 
the corresponding largest allowable tolerance within S, until we have 
entered r T , mt and its companion in the table. 
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An example should clarify these ideas. Ra for the voltage divider of 
Fig. la is reproduced in Fig. 3. Since the voltage divider has only two 
parameters, only one performance contour exists and this is identical 
to Ra- Suppose that Ri and R s have the same set of obtainable toler- 
ances, 

S t = S 2 = {15, 10,5,3, 1}. 

Figure 3 indicates L x ^ — 11 percent L 2 & 13 percent. We will there- 
fore consider pi to be the left side of the table. To begin to fill the table 
we notice m = 15 percent > |Li| and is therefore eliminated from con- 
sideration. t 12 = 10 percent is entered as the first entry in the left-hand 
column. Starting at the points in Fig. 3 (10 percent, 0) and ( — 10 per- 
cent, 0) , the vertical directions are explored to determine how far R 2 
may be varied before specifications are exceeded. In the second quad- 
rant, specifications are exceeded before A_R 2 = 2 percent. Therefore, 
t2,b = 1 percent is made the first right-hand entry in the table. Com- 
putationally we alternate between ARi = 10 percent and ARi = — 10 
percent and take small steps in &R 2 both positively and negatively un- 
til the specifications are violated (shown as step 1 on Fig. 3). The 
horizontal lines of the (10, 1) rectangle are then explored until Ai?i has 
been reduced to n,a = 5 percent. Since no specification failures are 
detected, the entry (10,1) is certified (step 2 in Fig. 3). (Notice we are 
assuming that performance contours are simply connected.) We then 
explore vertically from the four points on which step 2 ended and find 
that the next entry in the table is (5, 5) . We continue in this way until 



20% -t^ 2 (PERCENT DEVIATION) 




-20'/, 



1 *■ AR, (PERCENT 

20% DEVIATION) 



Pi 


P 2 


10% 


1%, 


5% 


5% 


3% 


5% 


1% 


10% 



+ -20% 

Fig. 3 — Allowable tolerance rectangles for voltage divider. 
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the table shown in Fig. 3 is complete, and we have explored the com- 
bined perimeter of all allowable tolerance rectangles. Notice that 
tightening A#i from 5 percent to 3 percent resulted in no increased 
allowable tolerance for -R 2 . 

This process is repeated for each of the Q) parameter pairs. The 
table for each parameter pair (r, s) is denoted T r ,. . Each table gives 
us a pairwise constraint on tolerances. The next step in our effort to 
find t* is to find the cheapest network satisfying the pairwise con- 
straints, i.e., to find a tolerance vector, t*, whose elements satisfy 
the table constraints pairwise and is the cheapest of all such tolerance 
vectors. Having found t*, we will then try to discover if r' t a r , 
i.e. if t* = t*. If t* i a, , we will find the next cheapest tolerance 
vector satisfying the pairwise constraints, etc., until we find one that 
is a member of a T . If we define a set, 

(jj T = { t | t satisfies pairwise constraints } , 

then it is clear that (B r 2 ft, and by choosing elements of (B r in order 
of increasing cost, the first chosen element of (B r that is also an element 
of & T is in fact r*. 

4.2 Branch and Bound Technique for Finding the Cheapest Network 

In this section we describe a variation on the branch and bound 
technique that is used to find the cheapest element of (B T , r*. Con- 
sider the multi-root tree structure shown in Fig. 4. The nodes in the 
first column represent the elements of & , the obtainable tolerances 
for p, . From each node in the first column we have branches leading 
to all the elements of S 2 , the obtainable tolerances for p 2 , etc. Finally, 
the last column has all the elements of S n repeated for every choice of 
t h e (n - 1) preceding element tolerances. There are II" -i m t nodes 
in the last column and an identical number of paths from the left 
side of the tree to the right. Each path is associated with a choice of 
tolerance vector and hence has an associated cost, C = 2^_, C^t*). 
The task of finding t* is then the task of finding a path through the 
tree from left to right such that: 

(i) The tolerance pair associated with every node pair of the path is 

consistent with the tabled constraints. 
(m) C is minimized. 
With each node we can associate a partial cost. For a node selected in 
column k, we have a partial cost C = £*-i C,(r,) where r k is the 
tolerance selected in column k and t, ,i < k, are the tolerances selected 
for the initial part of the path. 
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A flowchart for selecting the cheapest path whose associated rt(B r 
is shown in Fig. 5. In words, we begin at the left and choose for r x the 
largest (cheapest) tolerance not greater than Lj . We then move across 
the tree picking the fcth element's tolerance using 

r k = Min7\, t (r.), 
«<* 

where the function T iik ( Ti ) picks the right-hand entry opposite r f in 
the table, T iik . After each such selection we compute a partial cost, 
C k , for the portion of the path chosen so far. Two occurrences will allow 
us to eliminate that part of the tree with initial path n ,• • •, T fc _i , viz.: 

(t) At column k, one or more Ti, k {r t ) does not exist. This indicates 
that any t with the first (k — 1) components, t x , • • • , i-*., 
is not a member of (B T . 
(it) The partial cost C k ^ C T , a target cost. The target cost, initially 
set at co, is replaced by the cost of the first complete path 
found, and is updated as cheaper tolerance vectors emerge from 
the process. 




T ',2 



T '.3 




Tn,i 






REPEATED 


T n, 2 


DOWN THE 




- COLUMN 


T n,3 


fn-0 




XT mi. TIMES 


T n. 4 


L=l 



Fig. 4 — Tolerance choice tree. 



NETWORK TOLERANCE ASSIGNMENT 



1235 



SET t, = MAX {t,j| t,j<|L,|} 

k = i , C T = 00 




_I_ 






k = k+ 1 












1 








r k = min T L i, (r L ) 
L<k 






c t = i: o,(rj 



k = k-1 



1 F T k = fK . 

SET: 




MO 



YES 




YIELD T # = (T,,-..,T n ) 
AND C T 



NO SOLUTION 



Fig. 5 — Algorithm to find t* . 

In either case, we return to the (fc - l)st column and choose the next 
tighter allowable tolerance. If we are already at the tightest allowable 
tolerance in the (/c - l)st column, we return to the (k - 2)nd column, 
etc. If neither event occurs, a path reaches the right side of the tree. 
The associated cost of this path then becomes the target cost. We then 
back up to the (n - l)st column and choose the next tighter tolerance, 

etc. 

The efficiency of the algorithm derives from the fact that the vast 
majority of paths are never explored to completion, but rather are 
eliminated by one of the two terminating rules. 
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The algorithm terminates when all paths have either been explored 
or terminated. The final target cost is the cost associated with t". 
Having found t* , we must now determine whether t* t Ct r . 

4.3 Yield Checks Using Worst-Case and Tolerance Analysis 

A circuit built with the tolerance vector t* determined by the branch 
and bound optimization has the property that if any two parameters 
are allowed to vary within their tolerance limits while all other param- 
eters are held constant at their nominal values, design specifications 
will be satisfied. In other words, this circuit is the cheapest realization 
that meets all two-at-a-time constraints. Having determined that 
t* t (S> T , we must now verify that t* e CL T , or that the circuit with 
tolerance vector t* has 100 percent a priori yield. 

If it happens that t* i d T , we will return to the branch and bound 
algorithm, delete t* from the tree, and find the next cheapest network. 
We can thus view the branch and bound as a technique for ordering 
the elements of (B r by cost. At each call to the algorithm we obtain the 
next element in the nondecreasing cost sequence, t", r* 2 , • • • , t*°, • • • . 
The first element of this sequence that is an element of & T is r*. 

Tolerance vectors that satisfy the two-at-a-time constraints but that 
have less than 100 percent yield are easy to visualize. Figure 6 is a 




Fig. 6 — Ra and Rp for a circuit that meets pairwise constraints but has less 
than 100 percent a priori yield. 
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three-dimensional perspective of a circuit that meets the two-at-a-time 
constraints but has less than 100 percent yield. Considering only the 
positive octant, we see that each of the three performance contours is 
a triangle. After inscribing rectangles in each triangle, we have a cubic 
R P in three dimensions. Now if R A is a tetrahedron, we see that the 
shaded plane of intersection of R A and R P in Fig. 6 separates the pos- 
sible circuits into acceptable and unacceptable classes, and despite the 
fact that the pairwise constraints are met in each of the planes that 
includes the origin, we have a tetrahedron of failing circuits extending 
out of the drawing. 

Since we have not characterized R A , we go through two approxi- 
mate steps to determine if r* tft r . The first step is a quasi- worst-case 
computation. We use the signs of L, to extend each parameter to the 
limit defined by t, . In other words, we compute the circuit's per- 
formance with parameters 

Pi = P.o[l + T t sgn (L,)] 
where p i0 is the nominal value of p, . If this performance does not 
meet specifications, we know that the realization with tolerance vector 
t* cannot have 100 percent yield. If the worst-case circuit fails, we 
eliminate r* from the tree, go back to the branch and bound algorithm, 
and find the next cheapest circuit. 

Since even a linear circuit normally has a performance measure that 
is a nonlinear function of its parameters, the true worst-case per- 
formance need not occur at an extreme point of R P . Hence neither 
the preceding nor any other classical "worst-case" computation is 
sufficient to determine if R P C R A . The technique that has been 
chosen is to run Monte Carlo samples of the circuit using the tolerances 
of t* until either a sample fails, in which case we again return to branch 
and bound for the next cheapest network, or until 300 sample networks 
have passed specification. The probability distributions used for the 
Monte Carlo analysis are the triangular shapes shown in Fig. 7. This 
shape is chosen to emphasize the regions near the parameter extremes; 
satisfaction of the specifications when the parameters are near nominal 
has been assured by compliance with the pairwise constraints. 

After successful completion of 300 Monte Carlo samples, we claim 
that T* a t a T . We then replace the distributions of Fig. 7 with the 
actual distributions to be expected in manufacture and continue to 
run Monte Carlo samples as a further check. The latter process is 
continued to whatever extent is desired. It has been found in practice 
that some failing samples show up if the Monte Carlo analysis is run 
long enough. In this event we content ourselves with the notion that 
99+ percent is an acceptable yield. 
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V. AN EXAMPLE 

The techniques described above were tested on the bandpass filter 
whose schematic and requirements appear as Figs. 8a and b respec- 
tively. In our computation, 30 discrete frequencies were used to de- 
fine the specification. The network designer's first guess at a toler- 
ance assignment used all 2 percent components. In fact, such a choice 
does insure 100 percent yield, but it is not the optimum choice. 

The set of obtainable tolerances were given as: 

& = [20,15,10,5,3,2, 1,*}; i = 1, • • • , 8. 

Absolute costs, of course, are of no importance to the algorithm; all 
computations are performed on the basis of relative cost. Therefore, 
for this example we used 

C^) = 1/r,; i= 1, ••• ,8. 

It is interesting to consider the dimensionality of this example. We 
have eight parameters, each with eight obtainable tolerances. Hence, 
we have 8" = 16,777,216 conceivable tolerance vectors. Many of these 
vectors can be eliminated by not considering tolerances greater than 
the minimum intercept of the associated parameters. For example, 
the inductor labeled p 5 had L 5 = —5.7 percent. If this is done, the 
number of conceivable tolerance vectors falls to just under 3 X 10 8 
— still a rather formidable number. 

Two examples of the tolerance-pair tables are shown below: 



Ps 


P6 


5% 


Wo 


3% 


3% 


2% 


5% 


1% 


5% 


Wo 


5% 



Pi 


P* 


20% 


Wo 


15% 


5% 


10% 


10% 


5% 


15% 


3% 


15% 


2% 


15% 


1% 


20% 


Wo 


20% 



The cheapest network that the branch and bound algorithm located, 
i.e., that satisfied the pairwise constraints, was 

r* = {10,5,5, 10,3,3, 10,5}, 
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Fig. 7 — Probability density function used in Monte Carlo analysis. 

with a cost, C(t*) = 1.57. The efficiency of branch and bound is demon- 
strated by the fact that although there are almost 3 X 10 6 paths through 
the tree, r* is found in 200 ms of computer time on a CDC 3500. To 
understand this efficiency we note that the number of tolerance vectors 
satisfying the pairwise constraints is under 700,000. This is 4 percent 
of the original sixteen million realizations and less than a quarter of 
those remaining after each S { is diminished by {t { | t, > | L, | } . Perhaps 
more important to note is that the optimum realization has cost 1.5 



o — ^Wy\. 



P 3 



P., 



2 35 




240 360 490 700 

FREQUENCY IN HERTZ 



(b) 



Fig. 8a — Example bandpass filter. 
Fig. 8b — Example circuit specifications. 
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while any path that reaches a node of tolerance \ percent has a cost 
greater than 2. Clearly, a great number of paths are eliminated quickly 
because of the cost constraint. 

Examination of the algorithm's operation reveals that t* with cost 
1.57 failed the "worst-case" test. The branch and bound routine was 
then re-entered 2167 times to produce tolerance vectors of 25 different 
costs each of which failed the "worst-case" test. At a cost level of 2.47, 
362 tolerance vectors were produced of which three passed worst-case. 
However, each of these three failed during the TAP analysis. Costs 
of 2.50, 2.53, and 2.57 were then obtained with a total of 1124 vectors 
of which 12 passed worst-case but again failed during TAP. Finally, 
the tenth vector with a cost of 2.60 passed both worst-case, the initial 
300 TAP samples and a subsequent 500 TAP samples. The solution 
vector was: 

r* = {3,3,5,3,2,2,5,5}. 

The final cost of 2.60 compares favorably with the designer's first 
assignment, using 2 percent components, that has a cost of 4.0. 

The total time to execute the entire algorithm was 20 minutes on a 
CDC 3500. There were 8124 circuit evaluations: 121 to compute the 
L t , 5659 to form the tables T r ,, , and 2024 total TAP samples. 

VI. CONCLUSION 

In this paper we have described an algorithm that will assign net- 
work element tolerances in a manner that minimizes network cost. The 
tolerance for each element is chosen from a list of obtainable values, 
which, together with associated costs, is supplied by the designer. The 
set of tolerances chosen by the algorithm satisfies two conditions, 
namely: 

(i) All sample networks built with components whose values are 
within the limits imposed by the tolerances meet circuit per- 
formance specifications. 
(ii) The sum of element costs associated with the chosen tolerances 
is a minimum. 

The number of circuit evaluations required by the algorithm is rela- 
tively large and can be expected to grow as the square of the number 
of elements. An example circuit with eight parameters required just 
over eight thousand circuit evaluations. The feasibility of the al- 
gorithm is thus tied to the efficiency of our circuit analysis techniques. 
Linear networks, nonlinear static networks, 3 and a limited class of 
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nonlinear dynamic circuits are presently analyzed in a few seconds 
or less and hence are acceptable candidates for tolerance prediction. 
On the other hand, it is not unusual for the time-domain analysis of 
a nonlinear circuit to require several minutes of large scale computer 
time. The application of the tolerance prediction algorithm to such 
circuits must await algorithmic or hardware advances that bring the 
cost of circuit analysis down by at least an order of magnitude. 

The efficiency of the techniques is also influenced by the nature 
of the Region of Acceptability, R A , the shape of which is determined 
by the equations describing the network and by the performance 
objectives. More specifically, the algorithm depends on the adequacy 
of the characterization of R A by the performance contours. For the 
example circuit we noted that there were approximately 700,000 
tolerance vectors satisfying the pairwise constraints. However, only 
3663 of these vectors had to be tested before one was found that was a 
member of Ct r , i.e., before a vector guaranteeing 100 percent a priori 
yield was found. We deduce then that the performance contours were, 
for our purposes, a good characterization of R A . It is possible to 
imagine an R A for which this is not the case. In such a situation, where 
the number of vectors that are not members of a T but which still 
satisfy the pairwise constraints is very large, the number of circuit 
evaluations for both "worst-case" and TAP computations can become 
large enough to render the entire technique unfeasible. Fortunately, 
our experience indicates that in a large class of circuits, element pairs, 
e.g., LC and RC products and quotients, are dominant influences on 
circuit response, and that for this class of circuits the performance 
contours are an acceptable characterization of R A . 

All of this is to say that application of the tolerance specification 
algorithm is not cheap, and moreover may become quite expensive in 
some pathological situations. Hence, the algorithm should only be 
applied to networks whose anticipated production rate is large. In 
other words, one must balance out expected savings resulting from 
optimum tolerance assignment against the computational cost of 
achieving such an assignment. In any case, the variation of component 
cost with tolerance selection is seldom an issue for limited production 
networks. 

One final point should be mentioned, namely, the designer's ability 
to obtain accurate component cost data is essential to the success of 
the process. The generation of such data has often proved more dif- 
ficult than it might seem. For example, choice of a component of given 
tolerance for a large production circuit may in itself cause a pricing 
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change for that component. Further, when one is dealing with adjust- 
able elements for which an adjustment accuracy is sought from the 
algorithm, one is faced with obtaining costs on a process before it is 
implemented. The solution of this problem, of course, lies in the close 
cooperation between design and production engineering. 

It should be clear that the algorithm depends strongly on the as- 
sumption that only tolerance choices that allow 100 percent yield are 
acceptable. This assumption eliminates consideration of component 
probability distributions and the necessity of performing multi-di- 
mensional integration. It allowed us to use the performance contours, 
"worst-case," and TAP computations to eliminate large numbers of 
tolerance vectors and bring the problem within manageable propor- 
tions. 

The assumption of statistical independence of the components is less 
vital. In the algorithm as presented above, it meant only that the Re- 
gions of Possibility in two dimensions were assumed to be rectangular. 
Techniques that allow component correlation are subjects of future 
work. 
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